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in 
o 
o 

(N 
> 

o 

(N 

(N 
> 

oo 

o 

O 

o 



Oh. 



X 



S. Shankaranarayanan* and Musongela Lubot 
HEP Group, International Centre for Theoretical Physics, 
Strada costiera 11, 34100 Trieste, Italy. 
(Dated: February 1, 2008) 

The possibility that the scale-invariant inflationary spectrum may be modified due to the hidden 
assumptions about the Planck scale physics — dubbed as trans-Planckian inflation — has received 
considerable attention. To mimic the possible trans-Planckian effects, among various models, mod- 
ified dispersion relations have been popular in the literature. In almost all the earlier analyzes, 
unlike the canonical scalar field driven inflation, the trans-Planckian effects are introduced to the 
scalar/tensor perturbation equations in an ad hoc manner - without calculating the stress-tensor of 
the cosmological perturbations from the covariant Lagrangian. In this work, we perform the gauge- 
invariant cosmological perturbations for the single-scalar field infiation with the Jacobson-Corley 
dispersion relation by computing the fluctuations of all the fields including the unit time-like vector 
field which defines a preferred rest frame. We show that: (i) The non-linear effects introduce cor- 
rections only to the perturbed energy density. The corrections to the energy density vanish in the 
super-Hubble scales, (ii) The scalar perturbations, in general, are not purely adiabatic. (iii) The 
equation of motion of the Mukhanov-Sasaki variable corresponding to the infiaton field is different 
than those presumed in the earlier analyzes, (iv) The tensor perturbation equation remains un- 
changed. We perform the classical analysis for the resultant system of equations and also compute 
the power-spectrum of the scalar perturbations in a particular limit. We discuss the implications of 
our results and compare with the earlier results. 

PACS numbers: 98.80.Cq, 04.62.-|-v, 98.70.Vc 



I. INTRODUCTION 

Cosmological inflation is currently considered to be 
the best paradigm for describing the early stages of the 
universe [1, 2]. Inflation leads to the existence of a 
causal mechanism for producing fluctuations on cosmo- 
logical scales (when measured today), which at the time 
of matter-radiation equality had a physical wavelength 
larger than the Hubble radius. Thus, it solves several 
conceptual problems of standard cosmology and leads to 
a predictive theory of the origin of cosmological fluctua- 
tions. 

During inflation, the physical wavelength correspond- 
ing to a flxed comoving scale decreases exponentially as 
time decreases whereas the Hubble radius is constant. 
Thus, as long as the period of inflation is sufficiently 
long, all scales of cosmological interest today originate 
inside the Hubble radius during inflation. Recently, it 
has been realized that if inflation lasts slightly longer 
than the minimal time (i. e. the time it needs to last 
in order to solve the horizon problem and to provide 
a causal generation mechanism for CMB fluctuations), 
then the corresponding physical wavelength of these fluc- 
tuations at the beginning of inflation will be smaller than 
the Planck length. This is commonly referred to as the 
trans-Planckian problem of infiation [3-5] . 

Naturally, considerable amount of attention has been 
devoted to examine the possibility of detecting trans- 



Planckian imprints on the CMB [6-28]. Broadly, there 
have been two approaches in the literature in order to 
study these effects. In the flrst approach, the specific na- 
ture of trans-Planckian physics is not presumed, but is 
rather described by the boundary conditions imposed on 
the mode at the cut-off scale. In the second approach, 
which is of interest in this work, one incorporates quan- 
tum gravitational effects by introducing the fundamental 
length scale into the standard field theory in a particular 
fashion. 

In almost all the analyses performed so-far in the lit- 
erature, the trans-Planckian effects are introduced into 
the scalar and tensor perturbation equations in an ad 
hoc manner. In the standard inflation, the scalar and 
tensor perturbation equations are derived, from the first 
principles, in a gauge- invariant manner [29, 30]. How- 
ever in the trans-Planckian inflationary scenario, to our 
knowledge, such a calculation has never been performed 
and the scalar/tensor power spectrum, with the trans- 
Planckian corrections, were obtained from the presumed 
perturbation equations. 

With the possibility of detecting the trans-Planckian 
signatures in the current and future CMB experiments^, 
it becomes imperative to obtain the scalar and tensor 
perturbation equation from the first principles. (For re- 
cent work on the trans-Planckian constraints from the 
CMB, see Refs. [31-35].) In this work, we perform the 



^ For the current status and developments of WMAP and 
*E-mail : shankiSictp. trieste . it PLANCK, see the following URLs: http://map.gsfc.nasa.gov/, 

t E-mail imusoOictp . trleste . it http: / /astro. estec.esa.nl/SA-general/Projects/Planck 
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gauge-invariant cosmological perturbation theory for the 
single-scalar field inflation with the trans-Planckian cor- 
rections. The model we shall consider in this work is a 
self-interacting scalar field in (3 + l)-dimensional space- 
time satisiying a linear wave equation with higher spatial 
derivative terms. The dispersion relation [oj = uj{k)\ thus 
differs at high wave-vectors from that of the ordinary 
wave equation. Such a model breaks the local Lorentz 
invariance explicitly while preserving the rotational and 
translational invariance. The particular dispersion rela- 
tion we shall study in detail is 

u:^{k) = \k\''+h,\k\\ (1) 

where fen is a dimensionfuU parameter, bn < im- 
plies subluminal group velocity, while 6ii > implies 
superluminal group velocity. The above dispersion re- 
lation is a subset of a general class of the form ui^ = 
|fc|^[l -|- g{\k\/ko)], where (/ is a function which vanishes 
as A;o oo, and ko is a constant which sets the scale 
for the deviation from Lorentz invariance. It has been 
suggested that these general modified dispersion rela- 
tion might arise in loop quantum gravity [36, 37], or 
more generally from an unspecified modification of the 
short distance structure of space-time (sec for example 
Refs. [38,39]). Possible observational consequences have 
also been studied. For an up-to-date review, see Refs. 
[40, 41]. 

Even though the above dispersion relation breaks the 
local Lorentz invariance explicitly it has been shown that 
it is possible to write a Lagrangian for the above field 
in a generally covariant manner, consistent with spatial 
translation and rotation invariance, by introducing a unit 
time-like Killing vector field which defines a particular 
direction [11, 42, 43] (see Sec. (Ill) for more details). 
In this work, we use such a framework in-order to obtain 
the scalar/tensor perturbation equations for such a model 
during inflation. 

Using the covariant Lagrangian used in Ref. [42], we 
obtain the perturbed stress-tensor for the scalar and ten- 
sor perturbations about the FRW background. We show 
that: (i) The non-linear effects introduce corrections to 
the perturbed energy density while the other components 
of the stress-tensor remains unchanged, (ii) The non- 
linear terms contributing to the stress-tensor are pro- 
portional to fc^. Hence in the super-Hubble scales the 
trans-Planckian contributions to the perturbed energy 
density, as expected, can be ignored, (iii) The spatial 
higher derivative terms appear only in the equation of the 
motion of the perturbed infiation field {S(fi) and not in the 
equation of motion of the scalar perturbations ($). (iv) 
Unlike the canonical scalar field inflation, the perturba- 
tions, in general, are not purely adiabatic. The entropic 
perturbations generated during the inflation, however, 
vanish at the super-Hubblc scales. The speed of propa- 
gation of the perturbations is a constant and is less than 
the speed of light, (v) The tensor perturbation equation 
remain unchanged indicating that the well-know consis- 
tency relation between the scalar and tensor ratio will 



also be broken in this model. 

We obtain the equation of motion of the Mukhanov- 
Sasaki variable corresponding to the inflaton field. We 
show that the equation of motion derived from the gauge- 
invariant perturbation theory is not same as those as- 
sumed in the earlier analyzes [3-5] . Later, we combine the 
system of differential equations into a single differential 
equation in $ and obtain the solutions for the power-law 
inflation in different regimes. We also obtain the spec- 
trum of scalar perturbations in a particular limit and 
compare with the earlier results. 

This paper is organized as follows: In the following 
section, the theory of cosmological perturbations for the 
canonical single scalar field inflation is discussed and es- 
sential steps leading to the perturbation equation are re- 
viewed. In Sec. (Ill), we discuss the general covariant 
formulation of the Lagrangian describing the scalar field 
with modified dispersion relation and derive the corre- 
sponding stress-tensor. In Sec. (IV), wc obtain the per- 
turbed stress-tensor for the scalar and tensor perturba- 
tions. In Sec. (V), we obtain the scalar perturbation 
equation and the equation of motion of the Mukhanov- 
Sasaki variable. In Sec. (VI) , wc perform the classical 
analysis and obtain the form of the scalar perturbations 
($) in the various regimes. In Sec. (VII), we solve the 
perturbation equations in a particular limit and obtain 
the power-spectrum of the perturbations. Our results are 
summarized and discussed in the last section. In Appen- 
dices (A, B) we derive the equations of motion of the 
fields in the FRW and perturbed FRW backgroimds. In 
Appendix (C), we obtain the equation of motion of the 
Bardeen potential in our model. 

Through out this paper, the metric signature we adopt 
is {+,-,-,-) [44], we set fi, = c = 1 and l/(87rG) = 
Mp^. The various physical quantities with the over-line 
refers to the values evaluated for the homogeneous and 
isotropic FRW background. A dot denotes derivative 
with respect to the cosmic time (i), a prime stands for a 
derivative with respect to conformal time (77) and ^ de- 
notes a derivative w.r.t space components. We follow the 
notation of Ref. [30] to provide easy comparison. 



II. GAUGE-INVARIANT PERTURBATION: 
CANONICAL SINGLE FIELD INFLATION 

In this section, wc obtain the scalar and tensor per- 
turbation equations for the canonical single scalar field 
inflation. In the following subsection, we discuss key 
properties of the perturbed FRW metric and the "gauge 
problem" of the scalar perturbations. In the subsequent 
subsections, we discuss the matter Lagrangian and pro- 
vide key steps in obtaining the scalar /tensor perturbation 
equations. 
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A. Perturbed FRW metric 

We consider perturbations about a spatially flat (3+1)- 
dimensional FRW line element 



dt' 



^{t)dy? =a^ri){dr]^ -dx") , (2) 

where t is the cosmic time, a{t) is the scale factor and 
T] = J [dt/a{t)] denotes the conformal time. H is the 
Hubble parameter given by = a/a while Ti. = a' /a is 
related to the Hubble parameter by the relation H = H a. 

At the linear level, for the canonical single scalar fleld 
inflation, the metric perturbations {Sg^^i,) can be catego- 
rized into two distinct types — scalar and tensor pertur- 
bations. Thus, the perturbed FRW line-element can be 
written as 



ds^ = a^{r{)[{l + 24))dri^ -2^^Bdx'dr) 
- [(1 - 2^P)5,, 



(3) 



2^^^JE + h^j]dx''dx^\ 



where the functions 0, B, tp and E represent the scalar 
sector whereas the tensor hij, satisfying h] = d^hij = 0, 
represent gravitational waves. Note that all these first- 
order perturbations are functions of (77, x). For conve- 
nience, we do not write the dependence explicitly. 

The tensor perturbations do not couple to the energy 
density (Sp) and pressure (Sp) inhomogeneities. However, 
the scalar perturbations couple to the energy density and 
pressure which lead to the growing inhomogeneities. At 
the linear level, the two types of perturbations decouple 
and can be treated separately. 

The scalar and tensor perturbations have four and two 
degrees of freedom respectively. In the case of tensor 
perturbations, the two degrees of freedom correspond 
to the two polarizations of the gravitational waves and 
hence are physical. The scalar perturbations suffer from 
the gauge problem. (For a detailed discussion, see Refs. 
[29, 30].) However, it is possible to construct two gauge 
invariant variables, which characterize the perturbations 
completely, from the metric variables alone i. c.. 



a 



^ = 'tl)-H{B-E'). (4) 



Physically, $ corresponds to the Newtonian gravitational 
potential and is commonly referred to as the Bardeen 
potential while ^ is related to the perturbations of the 
3-space. For the single canonical scalar field scenario, we 
have $ = ^. 

$ and ^ are related to the pressure and density pertur- 
bations of a generic perfect fluid via the perturbed Ein- 
stein's equations. The pressure perturbations, in general, 
can be split into adiabatic and entropic (non-adiabatic) 
parts, by writing 

6p = cl6p + p'T, (5) 

where = p'/p' is the adiabatic sound speed [45-47]. 
The non-adiabatic part is ^Pnad = pT, and 

^_5p dp 



p' 



-p' 



(6) 



The entropic perturbation F, defined in this way, is 
gauge-invariant, and represents the displacement be- 
tween hyper-surfaces of uniform pressure and uniform 
density. In the context of canonical single scalar field 
inflation, only adiabatic perturbations are present, i. e., 
6p = c^(5p (where = 1) or F = 0. However, as we shall 
see in the later sections, the trans-Planckian inflationary 
scenario introduces entropic perturbations as well. 



B. Canonical scalar Held 

The dominant matter component during inflation is a 
spatially homogeneous canonical scalar field Tp(r]) (infla- 
ton). The Lagrangian density for the canonical scalar 
field <f{r],x) propagating in a general curved background 
is given by 



(7) 



where V{(p) is the self-interacting scalar field potential. 
The equation of motion and the stress tensor of the scalar 
field {(fi) in the conformally flat FRW background (2) are 
given by 



2Hip' -V^(p + a^ V^{(p) 



0, 



Ti:^''^ = ^''v^^<fi- 



-^"cp^pip-v{ip) 



(8) 

(9) 



where V^p = {dV{(p)/d(p) and refers to the Laplacian 
in the flat space. 

Let us consider a small inhomogeneous quantum fluc- 
tuations on top of a homogeneous and isotropic classical 
background. For the scalar fleld, we have 



(10) 



where one assiimes that the perturbation 6ip is small. 
The perturbed scalar field {S(p) and the perturbed stress- 
tensor of the scalar field {STj^^'^^), like the other scalar- 
type perturbation fimctions (4>, B^ijj, E), suffer from the 
gauge problem. (For a detailed discussion, see Refs. [29, 
30].) Similar to Eq. (4), it is possible to define a gauge- 
invariant quantity for the perturbed scalar field and the 
perturbed stress- tensor, i. e.. 



S'fi'-'''^ =6<fi + 7p{B- E'), 
<5To°'^'^ = 6TS + TS'{B- E') , 
5Ti^ =5Ti +jf(B-E'), 

dTi^^"^ = 5Ti +(t§- {B - E'),i . 



(11) 
(12) 



Separating the homogeneous and perturbed part from 
Eq. (8), we have 



Tp +2Wp +a^V^ = 0, (13) 
S^(9i)" + 2H<5^(3*^' - (14) 

+V^^ a" (5^(5^) - 4^'$' + 2y^ a^ $ = . 
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Similarly, separating the homogeneous and perturbed 
part in the stress- tensor (9), we get 



TO 



2 \ a 

(,gi) crpOiv) 



(16) 



C. Scalar and Tensor Perturbation equations 

In the earlier subsections, we obtained the 
gauge invariant variables related to the scalar field 
[S^piai) (m) STI^iv)-^ and metric perturbations ($,*). In 
this subsection, we outline the essential steps leading to 
the scalar and tensor equations of motion. Even though, 
this is a standard result, and can be found in numerous 
review articles (see, for example, Refs. [29, 30]), we have 
given here the key steps for future reference. 

From Eq. (16), it is easy to see that the non-diagonal 
space-space components of the stress-tensor are absent. 
This leads to the condition that $ = "P. Thus, the Ein- 
stein's equations for the perturbed FRW metric (3) in 
terms of the gauge invariant quantities are: 



^ a''5Tf^^ (17a) 



2 Ml 



2 a'5Tf^'^ (17b) 
a'STl'"", (17c) 



+ 3W$' + {2H' + H^)^ = 



2 Ml 



2M 

Pi 

1 -y .~i(gi) 



where STj^'^^'"'' is given by Eq. (16). The three perturbed 
Einstein's equations can be combined to form a single 
differential equation in $: 



■ n^ $ = 0. (18) 



The system of perturbation equations (18, 17b, 14) is 
quite complex. In order to extract the physical content 
more transparent, these equations are expressed in terms 
of two new variables - Q and u - which are linearly re- 
lated to $ and 5<^(5'\ i. e., 



Q = a 



7$. 



(19) 



Q is a gauge-invariant (Mukhanov-Sasaki) [29, 30] vari- 
able whose equation of motion is homogeneous. This en- 
sures that one can quantize Q in the standard way using 
the Lagrangian associated with its equation of motion. 
At the early stages of inflation where the quantum ef- 
fects are important, the equation of motion of Q helps 



in quantizing the fields and fixing the initial conditions. 
However, at the end stages of inflation where the rele- 
vant modes have crossed the Hubble radius and behave 
classically, it is easier to analyze the equation of motion 
of u. 

The equation of motion of Q is derived as follows: Sub- 
stituting S(p in terms of Q in Eq. (14), and using the 
relations (18, 13, 17b), we get: 



where 



g"_ v^g Q = 0, 

z 



(20) 



(21) 



The equation of motion of u is obtained by substituting 
the transformation (19) in Eq. (18): 



u" - V^u - —u = , 



where 



H 

a 



'-Xn^-n') 



-1/2 



2L 



(22) 



(23) 



Q is related to another gauge-invariant quantity C(= 

— {Ti/p )Sp+ip) by the relation Q = 2cg z(. The quantity 
( is time-independent on scales larger than the Hubble 
radius and, more importantly, related to the large scale 
CMB anisotropics (via the Sachs- Wolfe effect) [2]. 
Decomposing Q into Fourier space, we have 



Ms = 0> 



(24) 



where k = |k| and p^ = —Qk = ~2cszQ. The above 
equation is similar to a tiinc^-indcpendent Schrodinger 
equation where the usual role of the radial coordinate 
is now played by the conformal time whose effective po- 
tential is J7g = z" / z [48, 49]. The scalar perturbation 
spectrum per logarithmic interval can then be written in 
terms of the modes Pg as 



[fc' Vsik)] = 



27r2 ; 



(25) 



and the expression on the right hand side is to be evalu- 
ated when the physical wavelength (k/a)^^ of the mode 
corresponding to the wavenumber k equals the Hubble 
radius H~^. 

Before proceeding to the next section, we obtain 
the tensor perturbation equation in the FRW back- 
ground. As we mentioned earlier, the tensor perturba- 
tions hij{r], x) do not couple to the energy density. These 
represent free gravitational waves and satisfy the equa- 
tion: 



)p^ = 0. 



(26) 
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where /.iT = a/ife. This equation is very similar to 
the corresponding equation (24) for scalar gravitational 
inhomogencities, except that in the effective potential 
({/^ = a" /a) the scale factor a(r?) is replaced by z{rj). 



III. MODIFIED DISPERSION RELATION 
LAGRANGIAN 

In this section, we briefly discuss the general covari- 
ant formulation describing a scalar field with modified 
dispersion relation and derive the corresponding stress- 
tensor. 

As discussed in the introduction, to keep the calcula- 
tions tractable, we will assume that the scalar field with 
the high frequency dispersion relation is of the form 



u;^ = \k\'' + b,,\k\ 



(27) 



where > 0. The above dispersion relation breaks the 
local Lorentz invariance explicitly while it preserves ro- 
tational and translational invariance. Even though, the 
modified dispersion relation breaks the local Lorentz in- 
variance explicitly, it was shown in Ref. [50] that a co- 
variant formulation of the corresponding theory can be 
carried out by introducing a unit time-like vector field 
which defines a preferred rest frame. 



A. Covariant Lagrangian 

The action for a scalar field with the modified disper- 
sion relation takes the form [11, 50] 



-g (£^ + C^^^ + 



(28) 



where is the standard Lagrangian of a minimally cou- 
pled scalar field given by Eq. (7). The last two terms 
— ^coT ^.nd jCu — contribute to the modified dispersion 
relation of the scalar field. C^^^ corresponds to the non- 
linear part of the dispersion relation while describes 
the dynamics of the vector field u^. The two corrective 
Lagrangians have the form 



where 





-&11 




(29a) 


r — 


-X{g>'-'u^u^ - 1) 




(29b) 


F = 


- % , 




(30) 


dV = 






(31) 








(32) 



The covariant derivative associated with the metric g^^, 
is Vjii while bn and di are arbitrary (dimensional) con- 
stants. The tensor gives the metric on a slice of fixed 
time while T)^ is proportional to the Laplacian operator 



on the same surface. The fact that m'' is a rmit time-like 
vector (u^u^ = 1) is enforced by the Lagrange multiplier 
A. In the above, have the dimensions of inverse 

mass square and mass square, respectively while is 
dimensionless. 

The equation of motion for ip and obtained by vary- 
ing the action (28), respectively, are 



VV^i^i + V,^ = 26„ [V^(oVm''V.m'') 



(33) 



2diV.F"'' - Xu" = -6ii [V {V'^<fiu''d,<fi) (34) 
- VV^c^m'^V- V.(V''^w'')X)V] ■ 

Prom the above equation, we get 

A = b^^u^ [V^ {V^<fu''d^<f) - WV^ifu^V^ip (35) 

Using the results of Appendix (A), it is easy to show 
that for the FRW backgroimd Tp satisfies Eq. (13) which 
is same as that of the canonical scalar field. The field 
equation for gives A = 0. 

Before we proceed to the computation of the stress- 
tensor, it is important to know the current astrophysical 
constraints on the parameters of the model [40]: The 
constraints of the parameter di comes from the big-bang 
nucleosynthesis [51] and the solar system tests of general 
relativity [52]. These give: 



1 

M2 7 



(36) 



The constraints on the parameter b\\ comes from the 

observations of highest energy cosmic rays [53]. Using 
effective field theory with higher-dimensional operators, 
resulting in the modified field theory with a dispersion 
relation, it was shown that for various standard model 
particles 



6ii < 5 X 10" 



B. Stress tensor 



(37) 



In this subsection, we obtain the stress-tensor for the 
scalar field defined in Eq. (28). Formally, the stress- 
tensor for a general Lagrangian containing up to first 
order derivative in the metric is given by (cf. Ref. [44], 
p. 272) 



dC 



-g 



dC 



-.97 



didpgf''') 



■ (38) 



For simplicity, we will separate the contributions from 
the three different Lagrangians defined in Eq. (28), i. e.. 



cor) 



(39) 
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The strcss-tcnsor T^^,^ corresponding to the canonical 
scalar field Lagrangian is given in Eq. (9). The stress- 
tensor corresponding to the Lagrangian can easily be 
obtained and is given by 

+A \()p.v{ge^v!'u'' - 1) - 2u^u^\ . (40) 

However, the stress-tensor corresponding to C^^^ is much 
more involved. For the sake of continuity we give below 
only the final result while the steps leading to the result 
is given in Appendix (A). Wc get, 



(41) 



where, 



E^. = 4[9,[ln(V^)]Cpa.^ + ap(Cp5^(/p) 

- A^-M^^ - 5^^] (42) 



+ u^uP ^uJ'^uP - g^^uPu'' ^^P.u^u''] . (43) 

Before proceeding with the evaluation of the perturbed 
stress-tensor wc would like to mention the following 
point: Using the results of Appendix (A), it is clear that 
T\ji} and Tjfi,"^^ vanishes in the unperturbed FRW back- 
ground. Hence, the equations determining the evolution 
of the scale factor, i. e.. 



37^2 = _^T0^^^ ; 2n+n' 

PI 



(44) 



remain the same as in the canonical scalar field inflation. 
It is also worth mentioning that the trans-Planckian cor- 
rections do not play any role on the expansion of the 
FRW background while, as wc will see in the next section, 
the trans-Planckian corrections affect the metric and in- 
flaton perturbations. 

In this work, we will focus on the power-law inflation, 
for which, the scale factor is given by 

a{t) = {aotn or 0(77) = , (45) 



where p > 1, /3 < -2 (/3 
ao is a constant. 



/3 = 



2p- 1 
T 



P 



-2 corresponds to de Sitter), 



and (-?7o) 



(46) 



The scalar field potential and other background field pa- 
rameters are given by {q 



V = VoM^^ exp -q 



cp = aoM^^ In 



/2/p) 

; n -. 



ao 



V(/3 + !)(! + 2/3) 



z 
z 



i-v) ' 

V2/3(/3 + i) ; 

a/^ ^ PiP + 1) 
a (-77)2 



(47) 



IV. PERTURBED STRESS TENSOR 

In this section, we will obtain the perturbed stress- 
tensor for the scalar field with modified dispersion rela- 
tion (28) in the perturbed FRW background (3). As in 
the previous section, we will separate the contributions 
to the perturbed stress-tensor from the three different 
Lagrangians defined in Eq. (28), i. e.^. 



(49) 



The first term in the RHS of the above expression cor- 
responds to the perturbed stress-tensor of the canonical 
scalar field Lagrangian and is given by Eqs. (16). In the 
rest of the section, we will obtain contributions from the 
other two terms. 

Perturbing Eq. (40) and using the fact that Ff^,, van- 
ishes for the FRW background [See Appendix (A)], we 
get 

STi;(-^ = -2 6^6°{6X), (50) 

where SX is given by Eq. (B14). Using the fact that 
'D'^{ip) and dp'D'^{(f) vanish for the FRW background [See 
Appendix (A)], the perturbation of Eq. (41) takes the 
following simple form: 

JT^r' = 46n <5(DV) + C^^' OpiSV^^)] . (51) 

Substituting the relation (Bll) for d{V'^ip) and Eqs. (B7) 
in the above expression, we get 

- {rp" + mi) v^{5^)+iw^{5ip)''\ 51^51 , 

where ^ is defined in Eq. (B22). Substituting Eqs. (16, 
50, 52) in Eq. (49), we obtain the perturbed gauge- 
invariant stress-tensor to be: 



5T 



+ ■ 



4± 



1. 



\7'^^(9i) 



5T] 



iai) 



ST' 



o(s») 



(53a) 

S] , (53b) 

(53c) 



Following points are worth-noting regarding the above 
result: Firstly, the two corrective Lagrangians - £„ and 



^ The mixed stress-tensor 5Tlf is given by 

= S [g^'T^u) = gf^(STe.) + (Sgt^'yT:: (48) 
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C^^^ - contributes only to the perturbed energy den- 
sity {5p) . This imphcs that the non-diagonal space-space 
components of the stress-tensor are absent leading to the 
condition that $ = This also implies that the con- 
straint equation (17b) remains unchanged even for trans- 
Planckian inflation. Secondly, since the trans- Planckian 
corrections do not change the pressure perturbations, 
the perturbation equations for the tensor modes do not 
change. Hence, the tensor perturbation equations remain 
unchanged. Recently, Lim [43] had show that general 
Lorentz violating models (with out taking into account 
the higher derivatives of the scalar field) can modify the 
pressure perturbations and hence the tensor perturba- 
tion equations. However, in our specific Lorentz violating 
model, this is not the case. This indicates that the well- 
know consistency relation between the scalar and tensor 
ratio will also be broken in this model [60-62]. 

Thirdly, it is interesting to note that the trans- 
Planckian contributions to the energy density go as fc^. 
Hence as one would expect, in the super-Hubblc scales, 
only the canonical scalar field contributes significantly in 
these scales. Lastly, these can have two significant im- 
plications on the perturbation spectrum: (a) The speed 
of propagation of the perturbations (c^) can be differ- 
ent from that of the standard single-scalar field inflation. 
In the case of single scalar field inflation, we know that 
Cg = 1. However, due to the extra contributions to the 
energy density, this can no longer be true, (b) The per- 
turbations need not be purely adiabatic. ^ can act as an 
extra scalar field during the inflation and hence can act 
as a source. This can introduce non-adiabatic (entropic) 
perturbations. (See, for example, Ref. [54, 55]) We will 
discuss more on these in the following sections. 



V. SCALAR PERTURBATION EQUATION 



Perturbing the field equations (33, 34), we get 

5'^}'^''^ " + 2n 5<fi^'"'> ' - (S^^'"'>) +V,^^a'' 5^^"'^ (56) 

0, 



d„ 



Ml a2 



2di 



2M2j a2 



= 0, 



(57) 



where Ci = {M^^bii)/di is a dimcnsionlcss constant. 

Following points are interesting to note regarding the 
above results: (i) The spatial higher derivatives appear 
only in the equation of motion of 6ip and not in metric 
perturbation equation i>. Unlike the standard inflation, 
the perturbations arc not purely adiabatic and the speed 
of propagation of the perturbations is less than unity i. e. 
Cg = 1 — 2rfi/Mpj. (ii) In the case of canonical scalar-field 
inflation, the two dynamical variables $ and Sip^^"^"^ are 
related by the constraint equation (17b). In our model, 
$ and S(p^^^^ are again related by the same constraint 
equation, however ^ is related to the other flelds via the 
equations of motion. Hence, unlike the standard infla- 
tion, we have two sets of independent variables, (iii) In 
the case of canonical scalar-field inflation, $ and 
can be combined into a single variable — Mukhanov- 
Sasaki (Q) variable — in terms of which we can obtain 
the perturbation spectrum. However, in this model, as in 
the case of multi-field inflation models [54] , the equations 
of motion in terms of the Mukhanov-Sasaki variables are 
coupled. In the rest of the section, we derive the equation 
of motion of the Mukhanov-Sasaki variables correspond- 
ing to (5(^(3*^ and 

Substituting d(p in-terms of Q in Eq. (56), and using 
the relations (55, 13, 17b), we get 



Substituting Eqs. (53) in (17), the first-order per- 
turbed Einstein's equations take the following form. 



1 



2M2 
1 



2di 

t2 -, 



a 



(54a) 

(9i) 



(9i) 



(54b) 



L 

+ V,^aSip^^''>^ . (54c) 

As in the canonical scalar field inflation, the three per- 
turbed Einstein's (54) can be combined to give 



Q" 



where 



1 _ ) v^O -'-Q= ^V'Siv) , (58) 



2di 



Ml 



of 
n 



+ 



7,1/2 



(59) 



and Q^^', Q^^' are the gauge-invariant variables associ- 
ated with ^ and are given by: 



Substituting for ^ in-terms of Q^^^ in (57) , we get 



(60) 



+ 2(«-«|:)* = ||V««\ (55) 



3/2 



Q. 



(1) 



3-T_,' ^-T_,2 ^lA „2 
2^ - 2d^^ 



2Mla-'^ r n 



(61) 



, m = . 
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Decomposing Q, ,Q^^'' into Fourier space, we have 

_/2 "I 

fafc' 



^k'S.iv) (62) 



/if + 



2 4 2di a2 



n 



-3/2 



(63) 



ClH , 2 



2M2 a2 



where the Fourier transform of Q, Q^^\ Q^^\ respec- 



tively, are /^^ , /x^ , Q^^"* and 



(2) 



n 



(64) 



Eqs. (62, 63) are the main results of our paper, regard- 
ing which we would like to stress the following points: 
Firstly, in the earlier analyses, the equation of motion 
of the Mukhanov-Sasaki variable (Q) was assumed to be 
satisfy the differential equation Eq. (62) in which the 
source term was assumed to be zero. We have shown 
explicitly from the gauge invariant perturbation theory 
that, in general, this is not true. The RHS of (62) van- 
ishes in the super-Hubble scales (i.e A: ^ 0) where the 
perturbations can be treated classical. Hence, as ex- 
pected, the trans-Planckian effects are negligible. Sec- 
ondly, it is clear from Eq. (62) that the terms in the 
RHS will dominate during the trans-Planckian regime 
and can have interesting consequences on the primor- 
dial spectrum. Lastly, the perturbations (in general) are 
not purely adiabatic, i. e., it contains isocurvature per- 
turbations. However, these perturbations does not con- 
tribute significantly in the super-Hubble scales. Taking 
the Fourier transformation of the non-adiabatic part of 
the pressure perturbation ((5pnad), we have 

^(5Pnad)=-4rfl^Qf . (65) 

From the above expression, it is straight forward to 
see that, in the super-horizon scales, the entropic pertur- 
bations vanish. Following Refs. [45, 46], we can assume 
that, on large scales, the total curvature perturbation ( is 
conserved. As mentioned earlier, in the FRW background 
only the canonical scalar field contributes to the stress- 
tensor. Following Ref. [47], it is possible to show that, on 
large scales, only the curvature perturbation associated 
to S(j) contributes to the total curvature perturbation. 
Hence, it is sufficient to calculate the power-spectrum as- 
sociated to the scalar-field perturbation (Sip). This will 
be discussed in Sec. (VH) 



VI. CLASSICAL ANALYSIS 

In this section, we combine Eqs. (55,56,57) to obtain 

a single differential equation of $. We show that the re- 
sultant differential equation of $ is different from that of 



the standard canonical scalar field driven inflation. More 
importantly, the differential equation of <!> in our model 
is fourth order while in the standard canonical scalar 
field it is second order. We obtain the solutions of 
in three regimes — trans-Planckian (I), linear (H) and 
super-Hubble (III) — for the power-law inflation. In the 
following section, we obtain the power-spectrum of the 
scalar perturbations, in a particular limit, for the power- 
law inflation. 



A. The Power law inflation. 

Let us decompose the fields in their Fourier modes: 

mv,x) = aWe'^^-^ (66) 

We have dropped the superscript indicating that the 
quantities are gauge invariant. Combining the equations, 
we end up with a fourth order differential equation in the 
Bardeen potential. To keep the presentation light, the 

derivation of the equation is given in Appendix C. This 
equation, for the power-law inflation (45,48), reads 



(4) 



(67) 



+ 



•^2 \^ (_,,) (5+3/3) 



^1 + 



(_^)2{l+/3) 



S(2) 



re I 2 

(_^)3(2+/3) ^ T) ) (-r;)3+2/3' 



Tg Fio 2 
-T + -^k + 



(-'?) 



2,3+4 



(-'?) 



5+3/3 



The constants Fj depend on the background and the fun- 
damental constants in the following way 

Fi = 4(2 + /3), Fa = 12 + 13/3 + 3/3' - 6a2 , 
= ~2 ^1 rfT^^"' 

f_ 1)2/3 2+2/3 

r4 = 2 ^ ^ '° 6n, F5 = -18(l + /3)a^ 



Fe = 3 



al di 



M2, F7 = 6 + 4/3, 



/I \ 2/3 2+2/3 

Fg = -4(2 + /?) ^ 6n, 

= 6aX^2a!i_24 + ^^2^ Fio = 4 + 7/3 + 3/3^ 

C_l^> 2/3 2+2/3 

Fn = 2(l + /3)(2 + /3) ^ ^ &n , 



12 



3 (_ 1)3/3^3+3/3^2 



(68) 
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B. The zeroth order approximation 

We can find approximate solutions for the power law 
inflation in the following way. Introducing the quantity 
e by 



/3=-2-e, 



(e vanishes on the do Sitter space) we can make Taylor 
expansions of the coeSicients Fj and postulate the same 
for the Bardeen potential 



(70) 



m=0 



The outcome is the following. Each component ^k^miv) 
obeys a differential equation which is inhomogeneous, 
the source term depending on the preceding components 

As discussed in Rcf. [20], we have to be careful in the 
limit of e ^ in the sense that it does not give the per- 
turbation corresponding to the de Sitter space. This is 
related to the fact that on this background the inflaton 
field is constant so that the quantities Xi,Yi, Zi are un- 
determined. To work out what happens for the de Sitter 
space, one would have to consider the earlier equations 
where divisions by the derivatives of the infiaton was not 
performed yet. 

The zeroth order contribution is a solution of the equa- 
tion 



(71) 



We now introduce the dimensionless variable x defined 
hy X = krj and the function f{x) given by 



*fc,o(??) = fix) . 
The fourth order equation takes the form 



(72) 



r> (X) + ( 1 - A + fXx) - If'ix) + ^fix) = , 



(73) 



where, 



7 



= ^/2biiVoMpi , 



As mentioned earlier, we obtain approximate solutions to 
the above differential equation in three different regions. 



1. The first region 

In this region, the term 7^x^ dominates. In other 

words, the trans-Planckian effects are dominant and we 
are dealing with large values of x. Using the fact that 



f{x) = X is a solution of the full equation, one intro- 
duces the function h{x) by the relation 



f{x)=x J\{OdC 



(74) 



^59) and obtains the third order equation 



2(-l + 'y^x^)h{x) + 'y'^x^h'ix) + 4:h"{x) + xh" (x) = 

(75) 

We can get rid of the second derivative by the change of 
function 



h{x) = 



1 



R{x) 



(76) 



using the fact that we are in the region given by large 
values of x, one ends up with the differential equation 



R"'{x)+-f^x^R'{x) + WxR{x) 







(77) 



whose solution is a combination of generalized hyperge- 
ometric functions multiplied by polynomials: 



R{x) = C[Fp, |1 

-I- C'2s/^xFpq 

+ c's^x^Fp, 



1 3 

2' 4 



16 ' 



3 5 
4' 4 
5 3 
4'2 



^24 
J X 

16 ' 



16 



2 4 



where C-'s are constants to be determined and Fpg are 

the generalized Hypergeomctric functions. Thus, the so- 
lution to the differential equation (73) is 



fix) = Ciik)x 

+ C2ik)x^l^Fpq 

+ C3(fc)x^/3^P5 

+ Ciik)x^l^Fpq 







3 




1 


12' 6J 


l2' 


4' 


12 j 


' 16 



2 4 
7 X 



1 5 



1 



3 7 5\ _^ 2 4 

6' 12 J ' 1 4' 6' 4 J ' 16"^^ 



5 






17 




1 


12 




{1 


12' 


^} 


16 



2 4 

Yx 



(79) 



where C^'s are related to C-'s. These generalized hy- 
pergeomctric functions have a few properties which are 
worth mentioning. First, they are highly oscillating. For 
example, the function 





ri 3 




1 


12' 6 J 


l2'4' 


12 j 


16 



2 4 
J X 



(80) 



goes from 7.7 10*^^ to 1.4 10^®^ when x goes from x = 50 
to X = 100, fixing 7 = 1/10 for illustrative purposes. 

Let us now say a few words about these generalized 
hypergeomctric functions; this will help us to quantify 
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their oscillatory behavior. They are special cases of Mei- 
jer functions which can be defined by integrals on the 
complex plane [56]: 



where Ei{x) refers to the exponential integral. Using the 
asymptotic behavior of the exponential integral (cf. Ref. 
[57], p. 231), we get 



I ^1 h b2 ... bg 



where 



^^'^ = q_+,r(i-6,-s)n^_+ir(a, + ,s) ^^^^ 

and three possibilities are allowed for the contour 
C, according to some conditions on the parameters 
ai,bj,ni,n,p,q [56]. Our solutions correspond to m = 
n = 0. 

The asymptotic behavior which is relevant here is the 

following. For large values with — (i^* + 1)tt < argz < 0, 
the dominant part is roughly given by 



/-im,n I ^1 0,1 a2 ... \ rr (~„i-^tJ-*\ 

hi b2 ... bg J ^P"^^^ > ■ 

where 

H* = q — m — n, v* = —p + m + n, 

Hpq{z) = cx.p ({p - q)z~^z'' ,and 

q p 

I' 

q-p 



(83) 



p-q + 1 



In our case, one has to make the replacements 

z = -^l^kY, q = S,,p = 2 
lb 

so that 

3 



(84) 



(85) 



(86) 



—7 k ri j expl — 7 k ry exp(i7r/xj 



where cti = 2/3, (T2 = 5/3,(73 = 8/3. Fomr the above 
expressiom, it is easy to see that the solution in the 
the Bardeen potential $ is oscillating in this region. 
The choice of the constants Ci{k) correspond to differ- 
ent choices of initial conditions and thus, in principle, 
to different choices of vacua. We will come back to this 
later. 



+ Dsik) [-e-*^ + 2xsin(a;)] 
+ Di{k) [ie"-2xsin(a;)] (88) 

As we can see, the Bardeen potential is a sum of plane- 
waves. 



3. The third region 

When the term —Ijx^ dominates in the coefficient of 
the second derivative, the solution can be found and is 
given by 

f{x) = Gi{k) + G4{k)x + G3{k)x'^ + G2{k)xlnx{S9) 

From the above expression, we see that in the super- 
Hubble scales the scalar perturbations has a constant 
term which is identical to the canonical scalar field in- 
flation. 

To finish this section, let us remark that in the non 
trans-Planckian region, i.e when 



2 



1 - f2 » ^^^^ 



the solution to the difl^erential equation (73) can be ob- 
tained and is given by 



fix) = Hi{k)x + ^H2{k){2-2x + x^) 

1 / e~'^ 
-I- -H3{k)x[ iEi{-ix) 



1 / e*"^ 

- H/^ik) x I iEiiix) 



(90) 



this approximation covers the II and /// region simul- 
taneously. 



2. The second region 

In the intermediary region, 1 dominates over 7^0;^. The 
solution in this region is 

f{x) = Di{k) x + D2{k) x^ 

+ D:i{k) [e-'^'i-l + ix) -x'^Ei{-ix)] 

+ Di{k) [e"{i -x)+ ix'^Ei{ix)] (87) 



C. The first order approximation 

The first order contribution obeys the equation 



+ 2-2$fe,i(??) = 5fe(r?), 



(91) 
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where 



5 , 6„M>r 



pi' 

di 



10- 'r2,. ,42 , ,^2.. „2,_„/' ^7 



- —b,,M^iVok''r]^ + 4bnM^iVoV" log — 



+ 
+ 



24 5A;2 , „ 

6nM>r(-2di+M|) ■ 
+ 5^ ^ 



*fe,o(??) • 



This equation is exactly the one obeyed by the zeroth 
order contribution, except for the source term which is 
known since wc obtained the approximations of the ze- 
roth order in the three regions. Let us speciaUze to one 
of the regions and call Yi{ri) ,Y'2{rf) ,Ys{r}) ,Yi{ri) the four 
different solutions of the homogeneous equation given in 
Eq(71): 

The equation being linear and knowing the complete 
solution of the homogeneous equation ($ = '}2ti=i ^aXa 
with La constants) solution, we can solve it using the 
method of the variation of the constants. One can show 
that this can be achieved by the following system of equa- 
tions 



a=l a=l 



0=1 



0=1 



Let us concentrate on the second and third region for 
example (the non trans-Planckian zone) . One has 



Using the analysis discussed in this section, the power 
spectrum of the perturbations can be obtained upto a k 
dependent constant factor. In order to obtain the exact 
power spectrum, we need to quantize the theory and fix 
the initial state of the field [30] . In the following section, 
$^^^(ry)we obtain c;xac;t power spectrum of the perturbations in 
a particular limit. 



VII. POWER-SPECTRUM OF THE 
PERTURBATIONS - QUANTUM ANALYSIS 

In this section, we calculate the power-spectrum cor- 
responding to /Zg during the power-law inflation using 
the following approach: (i) Wc assume that the quan- 
tum field /Ltg is coupled to an external, classical source 
field «Sfc(77) which is determined by solving the coupled 
differential equations (55, 57). (ii) We solve the equa- 
tion of motion of ji^ in three regions - Trans-Planckian 
(I), linear (II) and super-Hubble (III) - separately [4]. 
We further assume that Sk{rj) will contribute signifi- 
cantly in the trans-Planckian region while it can be 
neglected in the linear and siiper-Hubblc region, (iii) 
The power-spectrum at the super-Hubble scales is deter- 
mined by performing the matching of the modes and its 
derivatives at the times of transition between regions I 
and II [(— ?7pi)^+'^ = (wfc^)"^/^] and regions II and III 
[r]jj = {1 + 0)/k]. We assume that the quantum field /Xg 
is in a minimum energy state sd r) = rji [58]. 

Region (I) corresponds to the limit where the non- 
linearities of the dispersion relation play a dominant role, 
i. e. 2bii[k / a{r])]'^ » 1 and kr] ^ 1. Region (II) corre- 
sponds to the limit where the non-linearities of the modes 
are negligible i. e. w ~ fc and krj ^ 1. Region (III) cor- 
responds to the limit where /cry <C 1. In the three regions, 
the equation of motion of /Xg (62) reduces to: 



(92) 



*fe,i(C) 



fcC / dr] 
/o 



(2i + 2^^--ifcV) . 
I bji(ikrj) 



+ 



{-2i + 2kr] + ik'^Tf) 



l-yfcC 
1 



2A;4?7 
1 



Ei{—ikr}) 



Skiv) 



-k^e 



1 



^ie~"'"^ + -xEi{-ikC) 



^ . (-2 + 2ikTi + li^rf ) 



drji 



lie^K - \ikCEiiikO^ (93) 

f' ^ A-2-2ik'n + k'^rf) , , 

X dri2'- ^ '"^Skiv)- 

A similar treatment can be applied to the trans- 
Planckian region but the formulas are too lengthy and 
will not be recorded here. 



t^f ^ -^k'S.iv), (94a) 



/3(/3 + l) 



(///) 



where 



l/2r_„„^(l+/3) 



'^'<^^ = '^0 5 ^0 = (25iJ / (-r/o) 



(94b) 
(94c) 

, (95) 



and Sk is given by Eq. (64) . The general solution to the 
differential equation (94a) is given by 

^^(t?) = Ar{k){-rtY'^Hi')[a{n)] 

+ A^ik) {-7iY'^Hi^)[a{'n)] + Mp(r?) (96) 

where Hpij]) is the particular solution to the inhomoge- 
neous part of the differential equation and is given by (cf. 
Ref. [59], p. 529) 
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ITV dlk^ 1/2 



Hi'\a{n)] [\- 

J ni 



s)i/'i/(''[a(s)]5fc(s)ds + 



Hl'\a{v)] [\- 



sy^^Hl^\a{s)]Sk{s)ds 



(97) 



r 



and 



2(3' 



a{ri) = ao{-v) 



ao = 



(98) 



^/ (< 111) is the epoch in which the integrals in (97) vanish. 
The quantities H^^ and H^'' in the above solution are 
the Hankel functions of the first and the second kind 
(of order v) , respectively, and the fc-dependent constants 
Ai(k) and A2{k) are to be fixed by the initial conditions 
for the modes at rji. Unlike the canonical scalar field 
inflation, where one assumes that the field is in a Bunch- 
Davies vacuum at r]i, it is not possible to assume such an 
initial condition due to the non-linearities of the modes. 
As mentioned earlier, we assume that the field is in the 
minimum energy vacuum state at rji, i. e., 



where /i^ (ry) = fjij^ (rj) — fip (rj) . From the above expres- 
sions it is evident that the particular solution to the in- 
homogeneous differential equation effectively changes the 
initial state of the field. Hence, flsivO can be treated as 
the new effective initial state of the field. It is worth men- 
tioning that the particular solution has not been fixed 
and is arbitrary. 



In Region II, the solution to the differential equation 
(94b) is (Minkowski) plane waves, i. e.. 



Ms(»7i) = 
We thus get 
Ai{k) = 

X 

A2{k) = 



; u (»?i) = ±^V 



VMm) 



(99) 



1 + 



(-ao/3) iVi) Hi%[a{r]i)] 

-^(-,o-/^A.(,oe\[.(..)] 

(-ao/3) i^AVi) Hl'}^[a{7j,l 



,(100a) 



(100b) 



)(j?) = Bi{k) exp[-ifcry] + B2{k) exp[ikv] , (101) 



where Bi{k), B2{k) are /c— dependent constants and are 
obtained by the junction conditions of the mode functions 
IJ,^p,li^J^^ and their derivatives at = rin. This gives. 



exp(— ifcT^pi) 

exp(ifc?7pi) 
(-r?Pi)V2 



Bi = 



B2 



ipap H^\[Q(r;pi)] 

iPao g^-'i[Q(r;pi)] 
fc(-?7Pi)'3+i jjm^airjPi)] 



1 + 



Hi'\a{r,Pi)] 



iBao gf_\[Q(?7Pi)] 
H-nPiy+' i/<''[a(,7p,)] 
ipao H^!ji[Q(??Pi)] 



1 + 



(2) 



a(??p 



, (102a) 
.(102b) 



In region III, the solution is 



M(iii)(??) = C{k)a{r])., 



(103) 



where C{k) is constant (not to be confused with the 
constants used in the previous section) whose modulus 
square gives the power spectrum of the density pertur- 
bations and is determined by performing the matching of 
the modes /x^"), /x^'") at r]„ = {1 + P)/k. We thus get, 



C{k) 



1+/3 



1+/3 



[Bi {k) exp{—ikr} jj ) + B2 (fc) exp(i/cr;^ )] . 

(104) 



The spectrum of the perturbations (25) reduce to 



(/? + l) 



47r2M2 p 



k^\C{k)\\ (105) 



We are interested in the leading order behavior of the pri- 
mordial power-spectrum and the possible modifications 
to the primordial spectrum due to the trans-Planckian 
effects. In order to do that, we need to obtain the lead- 
ing order behavior of the constants Ai,A2,Bi and i?2- 
Using the fact that krj ^ 1 and the asymptotic behavior 
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of the Hankel functions, viz. (cf. Ref. [57], p. 364) 



lim hI'/^\z) 



irz ' 



we get 



where 



Ai{k) « Aoils{Vi)eM-ixi){lTl) , (107a) 
A2{k) « ^oMs(%)exp(ia;i)(l± J) , (107b) 



1 - Mp(^i)/Ms('?») 

1 -Mp(%)/Ms(%) 



J = 



(108) 



Having obtained Ai,A2 in the Hmit of fcrj 1. Our next 

step is to cvahiatc Bi{k), B2{k) in the same hmit. In 
order to do that, we need to know the correct matching 
time 77pi. Demanding u)'^{r]p^) = k"^ gives (— Jjpi)^"*"^ = 
u)~'^/'^/k. We thus get. 



Bo 



nk 

-W 

irk 



1/2 



exp(ia;pi) 



1/2 



exp(— ia;pi) 



(109) 



(110) 



where ccpi = kr]pi{P + l)//3 — 7ri//2 — 7r/4. Thus, we get, 

2 



1 



Ms('7«) 



(111) 



X [1 + 2 cos(a;„ ) - 2Im[I\ sin(a;„ )] 



where 



Co 

X TT 



1 



{fi+l) 



167r2M2 /3 
2(1 + /3 - a;p, + Xi) 



1 + 



2(l+/3) 



(112) 



and we have neglected higher order terms Hke It is 
interesting to note that in the hmit of Skiri) 0, the 
power-spectrum is same as that of the standard power- 
law inflation spectrum with small oscillations. In this 
limit, wc recover the result of Refs. [4, 7]. In order to 
obtain the exact form of the power-spectrum, we need to 
evaluate jip which requires the knowledge of Sk{ri). 

In the rest of this section, we evaluate the power- 
spectrum in a particular hmit (1/ci 0). We, first, 
obtain the form of Sk{r]) by solving the system of cou- 
pled differential equations (55, 57) in two - sub-Hubble 
and super-Hubble - regimes. As mentioned earlier, the 
two differential equations (55, 57) do not contain higher 
order spatial derivatives. Hence, it is sufficient to obtain 
solutions in these two regimes. Performing the following 
transformations 



and taking the Fourier transform, Eqs. (55, 57), reduce 
to 



< + <k^ -^]^k = - 



2di e 



PI 

'2 fc2 



1 



Ci fc2 

a2 M2 



Ci fc2 

2a2 M2 



(114) 
(115) 



In the limit of 1/ci ^ (i. e. dilM^^ < b,,M^J, the 
above differential equations can be solved exactly. In this 
limit, the above differential equations become: 



^r=«(^«.)'- 



(116) 



For the sub-Hubble scales, during the power-law infla- 
tion, we get 

Ufe = Di{k) ex.p{—ikr]) + D2{k) exp{ikr]) 



2M2 



(_^,)(/3+i)Y/3(/3 + i) 
[Di{k) exp{ikr]) - D2(k) exp{-ikr])] (117) 

where wc have neglected the terms of the order 1 / (fcry) 
and Di{k), D2{k) are fc-dependent constants with the di- 
mensions of length squared (A;~2). Using the condition 
that the modes are outgoing, we set D2{k) = 0. In the 
super-Hubble scales, we have 



UkC^D3ik)a{rj) ; ^[^'^ = D^ik) 



2/3 2 



a'iv), (118) 



where Dz{k) is a constant. In the sub-Hubble scales, we 
have 



Sk{ri)=UMl\^ Di{k)k 



\ 13+1 

^ ) exp(-iA;7?) . (119) 
-V ) 



Our next task is to obtain [ip and the power-spectrum 

of the scalar perturbations. From Eq. (97) using the 
asymptotic limit of Hankel functions, we get 



£+1 , , 1-/3' 

^ (-r?o)^ 



X cos 



a(r/) -Fin ( r 



luMlk^l^ 



/3-1 
2/3 



(120) 



where we have set cx l/fc2. Substituting the above 
expression in Eq. (Ill), we get, 

{k''Vs(k)\ = Cofc'(^+')|l-Cifc(i+i/«|' (121) 
X [1 -I- 2cos(a;^) — 2/m[I] sin(a;^)] , 
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10 20 30 40 50 60 70 80 90 100 

FIG. 1: Plots of the standard power spectrum (thick curve) 
and the modified power spectrum with appropriate normal- 
ization. In plotting these spectra we have assumed that 
n = 10" GeV = 10^^ Mpc-\ CW/fcc) = 10"^ and p = -2.04 
where — k^^ (cut-off scale) and t^^ = Ti.~^ (inflation- 
ary energy scale). The above range of k/TC corresponds to 
2 < £ < 100 where £ denotes multipoles. 



where Ci depends on 6^^ and parameters of the power- 
law inflation. 

In Fig. (1), we have plotted the standard and the 
trans-Planckian inflationary power spectra. As can be 
seen the trans-Planckian power spectrum has oscillations. 
We would like to caution the readers that the oscillations 
are small, here we have magnified the efl'ect for illustra- 
tive purposes. We would also like to point the following: 
The power-spectrum (121) we have obtained becomes sig- 
nificantly different compared to that of Ref. [4] for very 
large k. For example, for a particular wave- vector this 
vanishes. However such an effect is not observable. 

Following points are to be noted regarding the above 
results: (i) We have obtained the general power-spectrum 
(111) of the scalar perturbations assuming that the scalar 
field is in the minimum energy state and that the contri- 
bution of the unit vector field to the energy density can 
be neglected. We have shown that the power spectrum 
depends on the form of the source term Sk which can be 
solved analytical in some particular limits, (ii) We have 
computed the power-spectrum of perturbations in a par- 
ticular limit i. e. (1/ci 0). In this limit, we recover 
the result of Refs. [4, 7]. 



VIII. DISCUSSION AND CONCLUSION 

In this work, we have computed the gauge-invariant 
cosmological perturbation for the single scalar field infla- 
tion with the trans-Planckian effects introduced via the 
Jacobson-Corley dispersion relation. Even though the 
dispersion relation breaks the local Lorentz invariance, a 



covariant formulation of the corresponding theory can be 
carried out by introducing a unit time-like vector field. 

Using the covariant Lagrangian, we have obtained the 
perturbed stress-tensor for the scalar and tensor pertur- 
bations around the FRW background. We have shown 
the following: (i) The non-linear effects introduce cor- 
rections to the perturbed energy density while the other 
components of the perturbed stress-tensor remains un- 
changed. Thus, for the trans-Planckian scenario, we have 
shown that $ = 5* and the constraint equation (17b) re- 
mains unchanged, (ii) The non-linear terms contributing 
to the stress-tensor are proportional to and hence in 
the super-Hubble scales, as expected, the contribution to 
the perturbed energy density can be ignored, (iii) The 
spatial higher derivative terms appear only in the equa- 
tion of the motion of the perturbed inflation field (Sip) 
while the speed of propagation of the perturbations [in 
the equation of motion of the scalar perturbations ($)] 
is different from that of the standard inflation, (iv) The 
speed of propagation of the perturbations (c^) is different 
from that of the canonical single scalar field infiation. (v) 
The perturbations arc not purely adiabatic. ^ act as an 
extra scalar field during inflation and hence can act as a 
source. This introduces non-adiabatic (entropic) pertur- 
bations, (vi) Since, the trans-Planckian corrections do 
not change the pressure perturbations, the perturbation 
equations for the tensor modes do not change. Hence, 
the tensor perturbation equation remain unchanged. Re- 
cently, Lim [43] had show that general Lorentz violating 
models (with out taking in-account the higher derivatives 
of the scalar field) can modify the pressure perturbations 
and hence the tensor perturbation equations. However, 
in this model, this is not the case. Since the tensor per- 
turbations remain the same, the well-know consistency 
relation between the scalar and tensor ratio will also be 
broken in this model [60]. 

We combined Eqs. (55,56,57) to obtain a single dif- 
ferential equation of <&. We have shown that the resul- 
tant differential equation of $ is different from that of 
the standard canonical scalar field driven inflation. More 
importantly, the differential equation of in our model is 
fourth order while in the standard canonical scalar field 
it is second order. We also obtained the solutions of $ in 
the three regimes for the power-law inflation. 

We have also obtained the equation of motion of the 
Mukhanov-Sasaki variable for the perturbed infiaton field 
with higher derivatives. In all the earlier analyzes, the 
Mukhanov-Sasaki variable (Q) was assumed to satisfy the 
differential equation Eq. (62) in which the source term 
(Sk) was assumed to be zero. More importantly, we had 
shown that the source term in Eq. (62) dominates dur- 
ing the trans-Planckian regime. The Mukhanov-Sasaki 
variable of the two fields are strongly coupled and hence 
obtaining the solution analytically is possible only in a 
particular limit. 

In this work, we calculated the power-spectrum corre- 
sponding to the infiaton field, during the power-law infia- 
tion by assuming that (i) the quantum field /ig is coupled 
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to an external, classical source field Sk (??) which is deter- 
mined by solving the coupled differential equations (55, 
57) (ii) the quantum field is initially in a minimum en- 
ergy state and (iii) di/M^^ <C b^^M^^. We have shown 
that in this particular limit, the power-spectrum is same 
as that obtained in Refs. [4, 7]. 

The work suggests various possible directions for fur- 
ther study: 

• We have obtained the power-spectrum analytically 

in the Hmit of rfi/Af^ < fe^^Af^. The trans- 
Planckian corrections in this limit are small to be 
observed in the present or the future CMB ex- 
periments. It would be interesting to obtain the 
power-spectrum by solving the system of differen- 
tial equations numerically and obtain the leading 
order trans-Planckian corrections in these models. 

• As we have mentioned earlier, this model intro- 
duces non-adiabatic perturbations which can lead 

to the non-Gaussianity in the CMB. Recently in 
Refs. [31-35], trans-Planckian constraints from the 
CMB was studied in detail. It would be interest- 
ing to do a similar analysis for this scenario. The 
non-Gaussian signatures may place stringent and 
independent constraints on the parameters b^^,di. 

• In this work, we have ignored the solenoidal part 
of the perturbed u field. The solenoidal part con- 
tributes to the vector perturbations. It would be 
interesting to see whether the solenoidal part of the 
perturbed unit-time like vector field can lead to the 
growing large-scale vorticity and hence the produc- 
tion of large scale primordial magnetic field. 

• In this work, we have ignored the back-reaction of 
the field excitations on the perturbed FRW back- 
ground. There have been claims in the literature 
[6, 9] that trans-Planckian modes may effect the 
evolution of cosmological fluctuations in the early 
stages of cosmological inflation in a non-trivial way. 
In Ref. [24], the authors have discussed in detail the 
backreaction problem of the trans-Planckian infla- 
tion in a toy model and have shown that the back- 
reaction of the trans-Planckian modes may lead to a 
renormalization of the cosmological constant driv- 
ing inflation.. It would be interesting to perform a 
similar analysis for this model. 

We hope to return to study some of these issues in the 
near future. 
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APPENDIX A: THE BACKGROUND. 

In this appendix, we give key steps in obtaining the 
stress-tensor corresponding to the two corrective La- 
grangians (29a, 29b), and the equations of motion of 
the scalar field and the unit vector field. Having ob- 
tained these, we discuss their properties in the FRW 
background. 

This question has been addressed in Ref. [11]. Our 
treatment differs with the one followed in that paper by 
the fact that we do not make the decomposition in time 
and space-like components. Our condensed formulas will 
prove very useful when computing the perturbations. 

In order to do that, it proves easier to go back to the 
action. Let us first specialize to the contribution of the 
non-linear part of the Lagrangian (29a): 

S^^^ = -K,j d'x^g{V^^f . (Al) 

Using the definition given in Eq.(31), the variation of 

V'^Lp can be written as 

+E''d^6>fi + H^^'d^dJ^ . (A2) 

The quantities A - ■ ■ H can be written explicitly. For ex- 
ample, 



Ajiu 



d 



HI"' = 



d 



DV- (A3) 



These partial derivatives have to be taken keeping in 
mind the choice of variables made in this work. In order 
to be consistent, we choose the following set of indepen- 
dent variables 

if, d^if, dadfjip, dadadfjip, gi"", dpg^"", drdpgi"', 

UcdfjUcdadfjUcX. (A4) 

Using the relations 



de{9.s9'')=^ ; ^ {9.s9'')=Q, (A5) 

09al3 



we get, 



Qgpa - "p"'-' Qgpa " 9po.9al3 ■ 
degaH = -9apg0<7deg'"^ ■ 



(A6) 
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The expressions giving the quantities A, - ■ ■ ,H are very 
long and time consuming to obtain and will not be dis- 
played here. One of the most important aspects of the 
way we will present our results is that for our purposes we 
only need to know their values on the FRW background. 
This is the subject of the next appendix. 
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Substituting Eq. (A2) in (Al) and integrating by parts 

the resultant expression, we obtain 



-fell / df^x 



-9{ -g^AT^'fY -2A^.V'ip 



+ < 



-26ii 
-26ii 



j (fx [-a^ (y^I^V^") (A7) 



From the above expression, it is easy to infer the 
contribution of the corrective Lagrangian to the stress- 
energy tensor as well as to the field equations of the in- 
flaton and the vector field u^, i. e., 



2&„ 



eq2,ip = — 26ii 



2 



(AlO) 



where (A9) gives the RHS of Eq. (33) and (AlO) gives 
the RHS of Eq. (34). 

Similarly, the variation of the action for the unit vector 
field (29b) leads to 



55, 



+ j d^x^{-g'^^u„U0 + 1) 8\ 

- drj d^x^gg^PgPP5{F^pFp,) . (All) 

From the above expression, we obtain the contribution 
to the stress-energy tensor and the field equation of the 
vector field, i. e.. 



Xu^u^ - 2dig'^''Fa^,Fp^ , 
1 



(A12) 



eq2,u = -2\g'"'u^-2d^^d,[^F'"'] , (A13) 

where (A13) gives the LHS of Eq. (34). 

It is easy to verify that the energy momentum tensor 

corresponding to the corrective Lagrangian and the vec- 
tor field vanish while the equation of motion of the vector 



field and the inflaton are satisfied when the following re- 
lation holds 



5^ = a(r/) (1,0,0,0) , A = , d^di"^ ■■ 



(A14) 



The intermediary relations needed are 

F^, = , T„/3 - - , (^) = . (A15) 



(A8) 



(A9) 



The evolution of the modified scalar field in the FRW 
background is same as that of the inflaton. Hence, the 
trans-Planckian effects, in this model, only affect the per- 
turbations. As we will see in the next appendix, the 
above equations play a significant role in discarding many 
terms from the perturbed stress-tensor. 



APPENDIX B: THE PERTURBED FRW SPACE 
TIME. 

In this appendix, we obtain the linear order perturbed 
stress-tensor corresponding to the two corrective La- 
grangians (29a, 29b), and the equations of motion of the 
perturbed scalar field and the unit vector field. 

We will use more than once the Leibniz rule: 

5{XY) = X5Y + Y5X = Xo5Y + Yo5X . (Bl) 

and the vanishing of the three dimensional Laplacian of 
the scalar field in the FRW background. 

For example, the variation of the component of the 
stress-energy tensor coming from the non-linear part of 
the Lagrangian reads 



^rpcorr 



g 



4A^^ 

- 4(s;,)^a,(5PV) • 



(B2) 



From the above expression, we see that it is enough to 
know the quantities A, - ■ ■ ,H on the background. Going 
back to Eq. (A2), we have to compute the variation of the 
vector field, the Christoffel symbol and other quantities. 
Let us begin by the connection coefficient: 



(B3) 



Due to our choice of independent variables, we need to 
express the variations of the covariant components of 
the metric in terms of the contravariant ones. This is 
achieved simply: 

5{gcc0g'^^) = Sgocn = -gajg/jaSg'^" . (B4) 
One then obtains 

<5r«. = aa' {-5-SX„+S^rj,.Sy5^rjpX 

+ 25^W.u - 25^<5°7?^<.) dgi'" (B5) 



pu 



+ 
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Similarly, we get 

+ (5oV + -5oV)^Wm- (B6) 

After a long but straightforward calculation, one obtains 
the rcsuh given in Eqs. (42,43 ) with the following ex- 
pressions for the background: 



In the same way, we get 

a" \ a a , 



(B13) 



{-^IJ,v)o 



2(f 



6 — (fi 
a 



(if)„ = -3(-7?'"^ + «) 



{E)o = 0. 



rip IT 



(B7) 



At this point it is worth noticing that — is the only 

non-zero component of the stress-tensor. Since the trans- 
Planckian corrections do not change the pressure pertur- 
bations, the perturbation equation for the tensor modes 
do not change. Hence, the tensor perturbation equa- 
tions is given by Eq. (26). Recently, Lim [43] had show 
that general Lorentz violating models (with out taking 
into account higher derivatives of the scalar field) can 
modify the pressure perturbations and hence the tensor 
perturbation equations. However, in our specific Lorentz 
violating model, this is not the case. 

Let us now specialize to the scalar perturbations. The 
covariant components read 



2(j) 



-diB 



-diB 2i)5'^ - 2didjE 



(B8) 



Using the identities given in Eq. (A5), we get the con- 
travariant components 



Using the fact that is a unit vector, we have 

5{g"^UaU0) = 0^5uQ = a(l) . (BIO) 

Using the form of the scalar perturbations Eq. (B8) and 
the results obtained in Eq.(B7), we obtain the following 
simple formula 

<5(DV) = ^ (^-^i^'s'^diSuj + V^S<fi^ . (Bll) 

Let us now turn our attention to the equation of motion 
obeyed by the perturbations. Concerning the inflaton 
field, the first two contributions to the formula given in 
Eq. (A9) vanish and the result is 



6eq 



26ii 



(B12) 



From the temporal index (/i = 0), we obtain the variation 
of the Lagrange multiplier, i. e., 

SX = 2di\5'^dj(Sui)' + ^\a<p'V'^5ip' -ip''^S'^diSu'. 

+ 5nip'^5'^djSui -a (ip" + 4W(^') V'^Stp .(B14) 
From the spatial indices fx = k, we obtain 

6ii (^-2a<p'\7'^dk5ip + 2{ip' f 5'^ dkdiduj'^ (B15) 



-I- Mio^ 



-Su'l + V'^Suk + dk{a(p)' - dkdiSui 



0. 



Having these results, especially the variation of the La- 
grange multiplier, we obtain the expression for the non- 
vanishing component of the stress-tensor: 

^j.(cor) ^ 2bn (^H^ng^s^^ _ ly2g^su', (B16) 



ST, 



00 



\ a a 
-2a2(5A. (B17) 



We now wish to express all the quantities in terms of 
gauge invariant quantities. For the vector field, we have 



(B18) 



5ui = 5u^i'''^ -uodi{B-E'). 
The equation of the vector field now becomes 

{a<p'V''dk5ip^<"^ - ip'^S'^dhdiSu'.'"^'^ (B19) 



while for the inflaton one obtains 

2V^^a''^ - 4^'$' + a^V^^^Sif'-'''^ + 2HS<p'^<''^' + S<p'^<''^" (B20) 



26ii 



0. 



Similarly, we can obtain expressions for the stress-energy 
tensor (53) interms of the gauge-invariant variables. 

Until this point we have not assumed any specific form 
of the spatial part of the perturbed unit-vector field. In 
general, the time-dependent spatial components of the 
perturbed u field can be expressed as a sum of irrotational 
and solenoidal parts, i. e. 



5ui = (V^i + (V X ro), . 



(B21) 
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The irrotational part of the perturbed u field will con- 
tribute to the scalar perturbations while the solenoidal 
part of the perturbed u field contributes to the vector 
perturbations. Since we are interested in the scalar per- 
turbations, for the rest of the calculations, we ignore the 
solenoidal part. Thus, we get 



Su) 



(B22) 



Substituting the above expressions in the field equation 
of the Bardeen potential, we get. 



$r+^3$fe'+(^2+i"2fc' + ^2fc^)< (C4) 

+ (Xi + Yik^ + Zik^)^k + {Xo + Yok^ + Zok^)<Pk = , 



Substituting this in the Einstein and the field equations, 
we obtain Eqs. (53). 



where 



APPENDIX C: THE FOURTH ORDER 
EQUATION 

In this appendix, we provide key steps in obtaining the 
equation of motion of Bardeen potential ($), in a general 
FRW background, by combining Eqs. (55,56,57). 

Prom the constraint equation we deduce that 



77% 



n 



(CI) 



Substituting this expression and its derivatives into the 
equation of motion of 5ipk, one obtains ^/j, i. e., 



where 



= 

Qi = 



3fell (^'2 F 



3&11 ip' 



if J ip 



1 

1 

fc2 ' 



+ 



36ii ^ 3 



36ii q) 

a 



n 



611 ip' 



1 

fc2 



36ii ip"- 



(C3) 
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